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ABSTRACT 


An approximate integral technique has been dev eloped 
for the hypersonic strong interaction flow over two-dimensional 
slender bodies with slip and temperature- jump effects. With 
the help of this approximate method, results have been 
obtained to understand the effects of slip and temperature 
jump on the various characteristics of flow such as heat 
transfer , skin friction, surface pressure distribution and 
boundary layer growth. The flow governing momentum and 
energy integral equations are solved by using a linear 
velocity profile and a second degree enthalpy profile, the 
induced pressure being described by the full tangent-wedge 
approximation. /. For both slip and no-slip cases, analytical 
solutions are dev eloped for the strong - interaction region 
by using the leading term of the tangent-wedge approximation. 
This is used to obtain the initial conditions for the 
numerical integration of the momentum and energy equations 
employing the full tang ent-wedg e approximation. The present 
results cover a wide range of flow situations with wall 
enthalpy ratios varying from 0.065 to 1, semi-wedge angle 
varying from 0° to 5°, Prandtl number varying from 0.75 to 
1 , and Mach number in the range of 10 to 20 for both slip 
and no-slip cases. Reasonable agr eement is observed between 
the present results and the experimental as well as theoretical 
predictions of other inv estigators . It is observed from the 
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present computations that except for the extreme cold wall 
case, slip and temperature- jump at the surface are very 
important in predicting the hyper sonic strong interaction 
flow close to the leading edge and in obtaining the physically 
plausible results . These results 3 further 3 help in bridging 
the gap between the experimental observations and theoretical 
predictions . 
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CHAPTER - 1 


INTRODUCTION 

§1.1. STATEMENT OP THE PROBLEM 

The problem of hypersonic rarefied flow past slender 
bodies has been the subject of extensive theoretical as well 
as practical investigations for quite a long time . In addition 
to its practical importance in the simulation of flight of the 
hypersonic vehicles in a rarefied atmosphere, such a flow model 
embodies a wide spectrum of flow regimes in gas dynamics from 
the kinetic flow close to the leading edge to the conventional 
continuum flow towards the downstream portions of the vehicle. 
Since the flow undergoes several transitions from leading tip 
to downstream, one will have to resort to different analytical 
methods to deal with each regime. 

At hypersonic speeds, the so called shockwave boundary 
layer interaction occurs when the rapid growth of boundary- 
layer displacement thickness near the leading edge displaces 
the external flow considerably. This phenomenon is the result 
of the low density of air due to the extremely high temperatures 
in the boundary layer at high Mach numbers. The displacement 
of the external flow gives rise to an induced pressure and 
imposes a negative pressure gradient on the boundary layer.. 
Simultaneously the shockwave emanating from the leading edge 
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must increase in strength to account for the added flow 
deflection, so that the stream lines crossing it carry 
increased values of vorticity and entropy. On blunt bodies 
the pressure induced by viscous interaction is generally 
insignificant in comparison to the bluntness- induced effects 
namely, the vorticity interaction. However, on sharp nosed 
slender bodies the viscous induced pressures may amount to 
several times the inviscid pressure. Also at hypersonic 
speeds the fluid density within the compressible boundary 
layer is relatively low which gives rise to the possibility 
of slip flow over the surface , 

A pictorial representation of the flow-field being 
discussed is given in Fig. 1, with the various regimes 
delineated. In both strong and weak interaction regions the 
surface pressures are found to be linearly dependent on the 
strong-interaction parameter x defined as (C/Re ) 7 . 

But as one moves further upstream the shockwave and boundary 
layer merge into each other and the pressure curve shows a 
11 plateau" although the slope of the shockwave continues to 
increase towards the leading edge. This levelling off of the 
pressure and of parameters like skin friction, heat transfer etc, 
are attributed to the shockwave getting thickened and weakened 
towards the leading edge. Theoretical investigations of the 
recent past could give explanations to these experimentally 
noticed plateaus when slip and temperature- jump at the surface 
are accounted for. The various diversified regimes up to and 
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including the merged layer regime occuring in the flow model 
have been treated by the continuum approach, the flow being 
governed by the Havier-St okes equations with the inclusion 
of velocity slip and t emp e nature- 3 ump boundary conditions at 
the body surface. Of course, in the region ahead of the 
merged layer regime the flow is highly rarefied and requires 
a kinetic theory treatment. 

§1.2. REVIEW OE LITERATURE 

Though the interest in the problem of hypersonic 
flow past slender bodies had been with the researchers since 
a long time > it was only in 1952 that a real contribution to 
the field occured when Shen £ 1^ attempted an analytical 
estimate of the viscous effects of hypersonic flow over an 
insulated wedge . Shen Q 2 suggested that the compressible 
boundary layer equations are valid in the viscous layer on the 
body as a first approximation and he was able to treat the 
shockwave boundary layer interaction problem with the help of 
a momentum integral method. Lees £3] solved the hypersonic 
strong interaction problem by using the Prandtl’s boundary 
layer equations. He found that as Mach number increases the 
interaction region spreads downstream over the plate and the 
skin friction coefficient increases rapidly over its convent- 
ional value. However Lees disagrees with Shen’s £ 2 ^ 
argument that the static pressure variation across the boundary 

layer is of the order (l/l^), for Shen puts up this argument 

1 1 

on the basis of the postulation that, p ^ ^ v —g , inspite of 
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the fact that I is certainly not of the order of near the 
outer edge of the viscous layer or near the surfaces for the 
case of heat transfer. Stewarts on £4~] extended Lees’ CO 
theory with a more complete discussion of the inviscid flow 
between the shockwave and boundary layer. He proposed the 
existence of a similar solution in the inviscid zone which 
may be joined to the similar solution in the boundary layer. 

He employed the von Karman-P ohlhausen integral method and 
showed that the tangent-wedge approximation for the pressure 
on the plate used by lees £ 3 3 is in deficit by ten percent. 

Pai and Shen £ 5 3 extended Li and Hagamatsu’s [[ 6 3 
non— adiabatic flat plate problem to v/edges "by using a generalized 
von Karman-P ohlhausen integral method. In their analysis they 
used a temperature function ’S’ to relate the absolute temperature 
and velocity. V/here, 


S = 1 - 


G p T + \ u2 

G 1 + ~ uf, 

P OO °° 


The major assumptions underlying their analysis were 


i) Pr = 1 
ii) u l n . 


Their results especially, those on surface pressure 
distribution were ’well supported by the similar solutions of 
Li and Magamatsu CSC- Betbram co extended Kendall’ s [I 8 3 
experimental investigations of hypersonic flow past flat plates 
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to thin wings with square and delta planforms. He found a 
greater value of the skin-friction coefficient which can not 
he explained by the conventional Jaminar boundary layer theory. 
Betram attributed this difference in skin -friction to the large 
values of induced pressures at hypersonic speeds. 

In the light of experimental investigations in a low 
density wind tunnel SGhaaf |~ 9 ^ suggested the existence of 
slip at the leading edge. In fact he was motivated by lees £lcQ 
who in one of his very substantial papers established the 
linear* variation of surface pressure with the strong interaction 
parameter. According to this result pressure must be infinite 
at the leading edge which is not expected to take place as 
the boundary layer assumptions are not valid in this region. 
Schaaf tried to modify lees’ theory so that the singularity 
in his formulation may be erased out. With the provision that 
the gas does not slip along the surface, such a modification 
would result in a surface pressure which is expected behind a 
normal shock. Therefore he concluded that the gas would have 
to slip along the surface invariably as the observed pressures 
were much less than what would have been behind a normal shock. 

Though quite a few experimental £ 11 ^ , [[ 12 []] , £ 13 ^ 
etc., as well as theoretical investigations £14^, £ 15^- etc., 
of the hypersonic slip effects have been made by many investi- 
gators, it was Talbot [J 16 ^ who first suggested the criterion 
for slip near the leading edge. Talbot concluded that the 
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parameter M (c/Re = V is appropriate for correlating the 

emergence of slip flow effects near the leading edge of a flat 
plate in hypersonic flow. He found that for V = 0.1 or above 
slip would be a dominant phenomenon. Talbot interpreted the 
rarefaction parameter in two ways. 

i) A measure of the ratio of the mean free path behind the 
shockwave to boundary layer displacement thickness. 

ii) A measure of the ratio of the leading edge shockwave 
thickness to distance downstream from leading edge. 

laurmann attempted a solution to the hypersonic 

slip flow problem over a flat plate by linearizing the Oseen’s 
equations, with an idea of getting some valuable approximations 
to an otherwise intractable problem as well as serving as a guide 
to methods of treating the full non-linear equations. His 
solutions showed that in the hypersonic limit the effect of 
slip and temperature -jump are quite radical, resulting in a 
flow at the leading edge which was unperturbed from free stream 
state. In otherwords there exists perfect slip at the leading 
edge. Aroesty £18^ with the help of a series solution argues 
that the slip effects alone are not adequate to explain the 
appearence of heat transfer plateau observed near the leading 
edge of a sharp flat plate in hypersonic flow as mentioned by 
Talbot. However, he agrees that the induced pressure and 
skin friction will exhibit a first order slip effect. Rudman 
and Rubin £19 3 developed a mathematical flow model incorporating 
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the effect of temperature- jump and slip at the surface. Their 
equations were valid throughout the boundary layer shock wave 
structure and inviscid oore for flat plates and wedges, for a 
wide range of Mach numbers >2. On the basis of the finite 
difference solutions, they put up the following conclusion^. 

i) Slip boundary conditions are necessary for V > 0.1. 

ii) Across the merged layer large gradients of pressure and 
density exists, the maximum density being on the surface'. 

iii) The merged layer passes asymptotically into the' strong 
interaction flow. 

Among the numerous experimental investigations, those 
performed by Vidal and Bartz • £2o^j are considered to be quite 
reliable and elaborate and therefore have been quoted in the 
literature extensively. They carried out a detailed study of 
the hypersonic low density flow over flat plates and wedges in 
a hypersonic shock tunnel at Mach numbers of 18 to 24 and 
Reynold’s numbers which produced conditions largely to span the 
transition between continuum flow to near free molecular flow, 
giving emphasis to determine the extend of slip and energy jump 
at the surface. His findings indicated that the then available 
theories substantially under- estimated the magnitude of slip 
velocity and energy jump. They discussed the viscous layer 
from a continuum view point and the first collision regimes 
from a kinetic view point. This model is an improved version 
of their own previous Innovation £21^. Their results are in 
good agreement with those of Shorenstein and Probstein £ 22 ^ 
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and Cheng [[ 13 3 for pressure and shock wave shape, hut heat 
transfer and skin friction are about 30 per cent larger than 
Cheng's. This error appears to stem from the approximation 
which Cheng has taken, namely, in the limiting case of V 
approaching unity 

P/Poo * M 2 sin 2 e, 
where e is the shock angle. 

In the vicinity of the leading edge,Yidal and Bartz found that 
free molecular limit is approached asymptotically with the first 
collision effects influencing the process. 

later, Kumar and Jain [[24 3 developed an approximate 
integral method for the hypersonic flow past insulated wedges 
incorporating the effect of slip velocity at the surface. The 
pressure along the edge of the boundary layer is prescribed by 
the leading term of the full tangent-wedge approximation and 
the von Karman-Pohlhausen method is used with a linear profile 
for the tangential component of velocity in the boundary layer. 
Chaudhury [[25 3 extended this work with full tangent-wedge 
approximation so that the solutions are valid from strong- 
interaction to weak-interaction zones. He solved the insulated 
we dge problem by feeding the attenuated values of pressure 
outside the boundary layer from the approximate solution. 

More recently Menon [[26 3 analysed the strong interaction 
problem incorporating the effects of slip and temperature- jump 
by a finite-difference technique suggested by Smith and Clutter. 
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Menon used an iterative procedure to obtain solutions to the 
slip flow problem from the no-slip series solutions of 
Ohattopadhyay Q27[]. 

§1.3. DESCRIPTION OP THE PRESENT METHOD 

In the present analysis an approximate integral method 
is developed for the hypersonic viscous slip flow over wedges, 
which is relatively simpler from computational point of view 
and gives reasonably accurate results. We envisage an infini- 
tesimally thin shockwave emanating from the leading edge and 
continuum hypothesis is assumed to be applicable to the flow 
field of interest namely the strong interaction regime, transi- 
tion regime and weak interaction regime. The results have 
been obtained by integrating the boundary layer equations (2.1), 
(2.2) and (2.3), using slip and temperature jump boundary 
conditions on the surface. 

Using the continuity equation a first integral of the 
momentum equation and energy equation is derived by assuming 
u * u^ and H Q v H^. Next, appropriate velocity and enthalpy 
profiles are to be chosen. The selection of the velocity and 
enthalpy profiles in any integral method is governed by the 
considerations of simplicity in the analysis on one hand and 
the possibility of satisfying more boundary conditions and 
compatibility on the other. A linear tangential velocity 
profile and a second degree enthalpy profile is found to be 
suitable for our purpose* The viscous- in vis cid interaction is 
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represented by the full tangent-wedge approximation so that 
the solutions are valid from the fore-regions of the body to 
downstream stations. The system of non-linear partial differe- 
ntial equations (2.1), (2,2) and (2. 3), is thus reduced to 
two second order ordinary differential equation (2 .34) and 
(2.43), which are coupled and solved by the scheme discussed 
in Ghapter III for the unknowns n , the boundary layer thickness 
and $ , the enthalpy gradient at the wall, lor no-slip case 
both these equations show singularity at the leading edge (i.e. 
at x = 0), hence the integration is started from a point slightly 
downstream. However, for slip case the integration is started 
right from the leading edge, lor both slip and no-slip cases 
analytical solutions are provided for the leading edge region 
by a series solution procedure employing the leading term of 
the full tangent-wedge approximation. 'This also provides the 
initial conditions for the numerical integration of the original 
problem employing the full tangent-wedge approximation. 
Integration is carried down to weak interaction region far 
downstream from the leading edge. 

Quite a wide range of flow situations have been 
considered with wall enthalpy ratios varying from 0 to 1, semi- 
wedge angles from 0° to 5°, Mach numbers from 10 to 20 and 
Prandtl number from 0.75 to 1 for both slip and no-slip cases. 

A detailed comparison of the present results with the available 
exper im ental data and results of other theoretical investigations 
have been carried out. For cold wall cases present results are 
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compared with the experimental findings of Vidal and Bartz [1 20 , 

Vidal and ffittliff £ 28 and Vas and McDougall £ 29 and a 
reasonable agreement is observed. Present results agree well 
with the locally similar solutions of Shorenstein and Probstein 
£ 22 3 and Cheng £ 13 and the finite difference results of 
Rudman and Rubins Q 19 ^ • ^ l ° r insulated v/edges a comparison 
is made with the approximate integral method of Kumar and 
Jain and for no-slip case present results agree well with 

Chattopadhyay ' s [[27^] results. Another comparison is made 
with the recent results of Menon £ 26 ] who solved the problem 
by iterative finite-difference technique. These comparisons 
show the effectiveness and accuracy of the approximate method. 
Similar to the earlier investigations, the present analysis 
also shows that slip plays an important role in explaining the 
experimentally observed 'Plateaus' of constant pressure, heat 
transfer, shin friction etc. towards the leading edge and 
bringing the result closer to the experimental data. 
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PCMiU RATION 01 TEE PROBLEM 

§ 2.1 G OVERNING BQ UATIOI S : 

The governing equations for a steady, two dimensional, 
compressible flow are the following : 

C^ontinu^it^ 


fe( p ") + 8j (pv) = 0 


( 2 . 1 ) 


M omentum 


„ 9u. 5u 

9x + M dx 


dP 


dP d__ 
dx + dy 


/ 5Un 

+ ~ y ay) 


dy 


= o 


En ergy 


n dH qH 5 

P U 4- PV ~r = ~ 

dx dy dy 


(2.2) 


JL ii + (1 _ i__) y (i u 2 ) 

Pr Sy + ^ Pr' y ^ ' 


dy ^2 


where H = 0 I + i u^ 
P 6 


B^uat^on of State 


P = P RT 


(2.3) 
(2.3a) 

(2.4) 


Also the viscosity-temperature relation p = v(T), for the 
present investigation is assumed to be linear. Accordingly, the 

following relation has been utilized : 

u T 

JE _ c „w 

y - 0 T 

■' CO OO ' ' 
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where c Is the Chapman-Rube s in constant. 

-§2 .2 BOUH DARY COIDIIIOHS 

(a) Boundary conditions for the momentum equation : 

At the surface of the wedge, for a no-slip flow 
situation, both the tangential and the normal components of 
velocity must be zero. However, for the strong- inter act ion 
region next to the merged-layer regime, the effects of slip 
will be important to account for. Accordingly, the boundary 
conditions have been given for slip and no-slip cases 
simultane ously . 

i) At the surface of the wedge, i.e. at y = 0 : 


u 

and 

u 

v 


v = 0, for no-slip case, 



C-, A 
1 w 


> 


v c)y'w 

for slip case 


where 



(2.6a) 


A is the mean free path of the molecules at wall 
w * 


temperature and is given by 


A = 1.256 
w 


t 1/2 


y /p a . 

w w w 


ii) At the edge of the boundary layer, i#e. at y = & : 


u = u_ 



14 


Based upon the hypersonic small disturbance theory, 
the tangential velocity at the edge of the boundary layer is 
equated with free-stream velocity, i.e. 


u = u e 1 Uoo (2.6b) 

(b) Boundary conditions for the energy equation : 

i) At the surface of the wedge, i.e. at y = 0 : 

H = (2.7a) 

where H-^ is the enthalpy at the edge of the Knudsen layer. 

ii) At the edge of the boundary layer, i.e. at y = 6 : 

H = H -w H (2.7b) 

The boundary condition (2.7a), in itself is incomplete, 
for the enthalpy at the edge of the Knudsen layer is still 
unknown. In the presence of slip it is observed that there 
exists a temperature jump across the Knudsen layer and the 
condition is 

h> - T w = °2 V (ff)w < 8 - 8 > 

where Gg is a. constant found to depend on prandtl number, 
the ratio of specific heats and accommodation coefficient. 

This dependence is expressed as 


> _ 2 Y iL 2-a 

'2 ~ y+1 Pr a “ 


(8.9) 


and, for a wall with full accommodation at the surface (i.e., 
a = 1) : 



15 


0 


2y 1 
2 “ Y+X p r 


11^ in (2.7a) can now be evaluated in terms of T^j (2.8) 
and u-^, (2 ,6a) . 

§2.3 SHE INTEGRAL EQUATIONS 

Using continuity equation (2.1), a first integral of 
the momentum equation (2.2) may be written as 

m i pu(u ~ - u) d y = H 6 + ( if)w (2 : 10) 


Similarly the integral energy equation is : 


dx 


/ pu(H-H) dy = - 


u. 


w / dH 


(§“X„ + (1 


Pr v by'w 


1_ X ,a__ 1 2x ' 
Pr ' v w'dy 2 u 'w 


( 2 . 11 ) 


§2.4 THE VELOCITY AND ENTHALPY PROFILES 


The distribution of the tangential velocity component 
across the boundary layer is expressed as 


u_ 

u 


= E 1 + k 2 y 


where k-^ and kg are constants. These constants, when 
determined by satisfying the boundary conditions at the surface, 
(2,6a) and the edge of the boundary layer, (2 ,6b) , yield the 


following linear velocity profile. 


u. = ° i x w + y # 

U °° °1 X w + 6 


(2.12) 
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The variation of enthalpy across the "boundary layer is represented 

by a second degree profile. This is compatible with the chosen 

2 

velocity profile as H v 0(u ). 


let , 


Ii 

H 


H = 7?“" = A + By + Gy" 


(2.13) 


where A, 3 and C are constants which are implicitly functions of x. 
Equation (2.13), when subjected to the boundary condition (2.7a) 
gives 

A ^ H 

1 = IT = H b 


Using the relations 


tt n m 12 

£ b = °p T b + £ u b 


H = 0 T 
w p w 


(2.14) 


and evaluating the gradients of enthalpy and velocity from the 
respective profiles (2,12) and (2,13), one can re-write the 
temperature- jump condition (2.8) in the form 


4 = s w - I a= x 


where, for a hypersonic flow 

2 


w 


u 


B 


0- X 
1 w 


“ H ~ (Gi 


(2.15) 




y 2 . 


Burther, as (~r) is a very small quantity, higher order 
terms have been left out. This approximation is justified even 


17 


for the extreme cases when H tends to zero, because a decrease 

w 

in would mean to increase B (= j ) and B becomes 

_ CO y 

appreciably large when H w reaches zero. Thus, in either cases 
(i.e. adiabatic as well as fully cold wall) the error introduced 
in the formulation due to this approximation is not appreciable. 


Thus. 


A = H + 0 o X B 
w 2 w 

By satisfying the second boundary condition (2.7b), 
one obtains 


G = 


7Z 


[ 1 - \ 


; - °2 X w B - 6 


Hence the enthalpy profile takes the final shape. 


H 


(Hw + °2 X w B ) + B y + 7S 

0 


H -B( i 
w v 


*2 A w + 6 ^ 


(2.16) 


How we have obtained the two expressions for the velocity 
and enthalpy profiles with 6 and B as unknowns. The coupled 
momentum and energy equations with substitution of these 
expressions are solved for these two unknowns. 

§2.5. SHE MOMENTUM EQUATI ON - A WORKING BORM 

Equation (2.3a) together with the equation of state 
(2.4) is used to express the density p in terms of P, H and u 
i.e . , 

TP 

(Y-l) (H - u^j 


P 


(2.17) 
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The momentum integral equation (2.10), after substituting 
for the density, p from (2.17) gives 


2 

d j 6 yP ^ u( 1 -u' 


dx J Y -1 H dy “ dx b + G 1 X " + b " 

Substituting for u and H from (2.12) and (2.16) and 
cancelling the common factor (1 - y/& ) from numerator and 
denomenator, the equation takes the form : 

2 


v u 

VV OO 


2YP6 


dx 




(°i x w + y) 


/ (- 

o (’-!)(% <V°2 x w ) + CV B (°a 


} dy 


dP 

dx 


u u 

W 00 


6 + ^7 A 


'1 A w n 

Phis equation is simplified further by noting that 

.2 




, nas~ « 1. 

2/ w%2 6 


(°1 X w +6) " (1 + -~ 2) 2 1 + 20 1 (^ 2 ) + C 2 (^) 

(2.17a) 

The above approximation has been made to facilitate the 
integration and we get, 

(°i \ + y) 


d r f 2YP 

{ 1 1^7 7s 


(V°2 x w B)6+CV B ( Q 2 V & >3y 

dP 


} dy 


dx 6 + Tcrv + b y * 

In this equation, P, B, and 6 are functions of x only 
and hence , these may be taken out of the integral sign. The 
resultant equation may now be integrated to yield 


v u 

w 00 
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2YP 


dx 


( r - 1 )CV a 2 V^B 


6 - 


V&-Ol V+(0l-0 2 )BS V 

\ + B (°2 X W + 6 ) 


3^ ( 2(5 w + C 2 \ B) + B6^ 


} 


dP . U q° 

dx 6 + (0 1 X w + 6") 


< H w + °2 X w B ) 

(2.18) 

The logarithmic term in the above equation is simplified 
by expanding it in a series form and neglecting the higher order 
terms. Thus, after some algebraic simplification the momentum 
equation (2.18) is simplified to : 


d 2YP6 
dx 1 Y~r” 


1 


S( V°8 X w 3 > +B6 


dP o v w u °° 
' } 6 + ?f“T 


1 w 


+ 6 


(2.19) 


Equation (2*19) is next differentiated and the coefficients of 
derivatives are grouped together, leaving terms of order 

(g ~f, g“ f“ and smaller. 


P { 3(i w +0g X w B)(g) - 6 2 (g) > + 6 Q 3(H w +Cp X w B) + B6 J 


'dx' 


. T H -l. 

w 2 w 


(1 -5T C3(l + 0„ > B) t 16]) <£) 


2 Y 


"w T w 2 w 


dP 

dx' 


Y— 1 

ST 


u u 
w °° 


L 3 ( \ + °2 X w B ) + 36 3 + 0 1 'TJ^+ & (2.20) 


A particular case of equation (2.20) for adiabatic wall 

- 2 

is obtained by putting H w = 1, B = 0, and by putting ^ Y % 1, 

one gets 
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■+* 6 


dp _ 
dx “ 



u 

CO 



( 2 . 21 ) 


In the method developed by Kumar and Jain j^24^ for 

the flow over an insulated wedge, the above form of the momentum 
equation was used. 


fhe range of validity of the equation (2.2o), depends on 
how the pressure and its derivatives are prescribed. We first 
use xhe leading term of the tangent-wedge approximation, which 
is valid for the strong interaction region only j 


M 2 (p + ||) 2 

or P = Poo u 2 (p + ||) 2 

and H = ( Y +l) P., u 2 (p + g) £-^ 

dx 


Substituting the expression contained in (2.22) intc 
equation (2..2o), one obtains, 


3(H +0p X b) & - & 2 P 

v \rJ A w / ni -v- y 


(f> + 


w 2 w^' as 

3 (\ + C 2 X w B) + B6 


3( V C 2 V) + BS1 


-I 


{ 1 Y “ 1 
11 - rr 


} ( B + ^ ^ 
i '•P + ’“■JJ 

^ dx 


y-1 

YfY+TJ 


5 ^w + G 2 X w B ^ + 136 


2 y 


P ~u orr + 5 (2.23) 

00 1 Y'l u 


Kext the various variables are non-dimensionalized as 


given below : 
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We first notice that the ratio, u / p u , has the 

7 W °° 00 7 

dimension of length. This ratio is denoted by 1 and has been 
used for non— dimensionalising the various quantities of interest. 


and 


£ = 


6 

I 

x 


d6 

dx 


4> = BL 
dn t 

= M = n 


2 2 

d!5 1 dn 

TTT * I 77 
dx d£ 


i ri n 

1 


dB 1 dv 1 ,,,> 

dx = d? = 5 ? * 


B6 = n<i 


(2.24) 


Also, let us define the two strong- inter act ion slip 
parameters and e^g as : 


'll 


G. X 
1 w 

T 


(P 


V ) 2 = 


/ 1.256 \ 


2(y-1) H 


w 


1 / 2 ' 


'e 


22 


(\— ) (P. 


» y 2 

n ) = 


/ 1.256 s p 
(" y 7 T “ '°2 




Y \ 


1/2 


(2.25) 


Equation (2.25) shows that under an extreme cold wall 
situation (i.e. 1 7 = 0 ), the phenomenon of slip is apparently 
absent. Insertion of the expressions (2.24) and(2.25)in eq.(2,23) 
leads to 

to. No. A 55AU * 

•• • 


f See appendix 



C3(H w +* 22 ) n - n Z *' U(p+rf)+nn" C 3 (H w + *gg) + nt H 


{2 ~¥ L3(\H.* g2 ) + n*]} 


Y— 1 


= yjf+TJ ZM\ + * 22 ) + x 


(P + r i’ ) 


whe re 




e 22 ♦ 


22 (P + n')^ 


E e n + n (P+ Ti' ) 2 |] 

(2.26) 

(2.27) 


Equation (2.26) gives a particular case of adiabatic 
wall with il = 1, and ip = o, like 


e nC '>'0'+|3)+nn'';] + nn' ( n’+p f + Z „" ( n'+fs) 2 = ^ ( „'+p) 


(2.28) 

This was the equation used by Kumar and Jain £24(] for 
the analysis of hypersonic flow over insulated wedges with e^ 
here identically equal to e. 

The pressure and' its derivatives In the equation (2.2o) 
are now represented by a full tangent-wedge approximation 
thereby, extending the range of validity cf the analysis to 
weak- interact ion zone also. The full tangent-wedge approximation 
is expressed as 



23 


or 


2 

p u 

CO CO 


p = ~“F~" { 1 + (vk) [IT 1+?1? + lok^} 


YM 


and 


p n‘ 

(3JP oo oo 


I = uo_ — , ( £|) 

OO n j 

T i + ff dx 


(2.29) 


where 


k = 


and b = 


M ~ « + S> 


Y+l 


”4” 


(2.30) 


Here also we define the following two slip parameters e^ and e { 


e. = 


1.256 G-^ ~\f ■ . T H 


and e g = 1.256 C 2 f § w 


(2.31) 


Also we can see 


and 


+ C, X 
1 w 

”“IT” = 

+ 0 *■ 

°2 w 

““IT” = 


e 1 w 


e 2 “ 


(2.32) 


where 




0) — 


P/F 


See appendix 
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or, from equation (2.29), 

“ = + yk f TT1? k^ + b k]> (2.33) 

Equation (2.29) through (2.33) are now substituted in 
equation (2.20) after non-dimens ionalising it with the help of 
relations given by (2.24). Thus we obtain final working form 
of the momentum equation which incorporates the effects of both 
slip and heat transfer. 

H 3 ( H w +e 2 t0 ^) I] f l+Yk£ T l+b 2 k^+bk "J } + nn" £ 3(H w +e g ijrio)+ rnjTJ 


{ 1_ { 2y mJ^Y 1+bV + bkj 


:t"z 


—r?-7? 

Y l+b k 


+"*3 2 7— - 1 - 


( n + e 1 to ) 


(2.34) 


A particular case of equation (2.34) for insulated 
wedges with H w = 1 and ip = 0 results in an equation which is 
exactly similar to Ohowdhary’s £25^] equation with full tangent- 
wedge approximation for adiabatic surfaces. 


§2.6 211 ENERG Y EQUATION - A W ORKING M 

The energy integral equation (2.11), after substituting 
for the density from (2.17) and the gradients of enthalpy and 
velocity from (2.12) and (2.16), takes the form 


d A YP 

di J Q TfZTj 




(2.35) 
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Substituting for u and H from (2.12) and (2.17) and 
neglecting terms of order (~£) 2 * * * or smaller, equation (2.35) 
is transformed after the algebraic simplifications to a form 
more amenable for a closed form integration in the normal 
direction as, 



*00 

TCT)6 


C \ + ( 0 2 \ + S ) B I] 


(Z. 

y 


-+-2E)} 
+ n 15 


dy 


where 


m 


and 


n 


.w r u 
P r B 


2 n \ 

1 u 0 1 A 
/ l \ 00 1 w 

“ Pr' 


(°i x w + &) 


( 5 W + c 2 X w B - 1 ) & 

Cy»(« 2 X w + 6)-l] 

+ C 8 V & 

Cl, + (Og \» + £■) - 


(2.36) 


(8.37) 


After carrying out the integration equation, (2 .36) 

yields : 




yp u. 


ZK + (°2 X w + 

+ (Og + 6) Bp 


6* 

F 


+(m-n)&+n(n-m) ln(^™) 


pf B 


2 ~ } 

. u G* A 

... 1 . . eo r 1 W 

— (^ " pr^ x rr^ - 

(°1 w + 6) 


(2.38) 


Ihe logarithmic term in the above equation is simplified 
by expanding it in a series form and neglecting the higher order 
smaller terms : 
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3ja ( 


8+n 
"6 ' 



{ 5L + ( °2 X w + 6) B 

3 (\ + C 2 X w B ) + 26 


(2.39) 


Upon insertion of equations (2.37) and (2.39) in eq.(2.38) 
and after considerable algebraic rearrangements and simplifications, 


We obtain for equation (2.38), 


d , YP u oo & 

3x { ^TTSfT 


3 (\ + °2 X w B ~ , 1 ) + 36 } = %r 

-3(C + c 2 v;bT = ' Er 


H 


- (1 - 


Pr 


) v u - 

w 03 (0 X X 


c. * 

1 w 


w 


+ 6 ) 


(2.40) 


Carrying out the different iat ion in equation (2.40) and grouping 
the coefficients of derivatives together, one obta,ins, 


sxvrrj { P L { C 3 (\ +G 2 x w B ) + 35 3 2 _ + °2 x w (£) 

+ 36 (3C 2 X w + 6) §: [3(V°2 6 


H 3(\ + °2 x w B ) + 36 3 } 


y w 

Pr 


- (1 - fr) 


. u 2 
00 


Ci x 

1 W 


(c 


1 w 


+ 6 y 


(2.41) 


The pressure and its derivatives appearing in eq. (2.41) 
are first described by a tangent— wedge approximation valid for 
the strong interaction region given by the expression (2.22). 

The equation is non-dimens ionalised by employing a scheme similar 
to (2.24). The slip parameters employed here are the same as 
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those given by expression (2.25). Thus with the substitution 
of equations (2.22), (2.24) and (2.25) in equation (2.41), we 
obtain the energy equation for the strong interaction region 
in its final form as : 

(p+n') 2 < n'[(3(^ + * 2g ) + n*} 2 - 9(1 + * g2 )3 + 3 n *' (fl— 

(P+ n ) 

+ 2 nn" (p+n') C3(V*82) +n *3 C®(V*82) +n *- 3 3 

♦ , e ll(P+"') 2 1 

%Pr “ Pr ' ~ , ' 

e 11 + n (p+ n ) -J 

(2.42) 


= C3(V*g2)+"€ 2 


where is defined as given by (2.27). 

We now extend equation (2.41) with the help of full 
tangent-wedge approximation given by expression (2,29). Non- 
dimens ionalising the equation (2.41) in the same way as 
equation (2.24) we obtain 


[(" 3(iI w +egO)ii» )+rnjQ ^ - 9(H y? +eg“ 1 i') 3 n +3n (3egW+ n) ^ 


{1+YhQT T+*??+bk3 } + nn" il 3(H w +e 2 # )+ *<T| C 3( V e 2 
+ nt -33 { 2Y M m cr l+bV+bk 3 - — ° — > 


= - 2(y-1)i 4L 3CVe 2 ^) + n<T 


~~ 

f i+b V 

e . ^ 

* +( i i)— n 


• E ’ r/ (e^+n)^. 


(2.43) 


+n )) 
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where k is given by (2.30), and e ± and e 2 are given by (2.31), 
and co is given by (2.33). 

Equation (2.43) represents the final working form of 
the energy equation which incorporates the effects of both 
slip and temperature-jump and is valid for the strong interaction 
region as well as for the weak interaction zone. ' 

Thus, we have the two equations, namely ^the momentum 
(2.34) and the energy (2.4*), v/hich are coupled and solved 
simultaneously for the unknowns n , the non-dimens ionalised 
boundary layer thickness and i|> the non-dimensionalised 
enthalpy gradient. 

§2.7. WO RKIN G- EQUATIONS NOR NO- SUP CAS E 

The governing equations for the no-slip case are 
derived as a particular case from those for the slip case, with 
the slip par am eters equated to zero. Thus the momentum and 
energy equations for strong interaction regions for no slip 
case are obtained from (2.26) and (2.42) by putting e^ = egg = 

*22 = °* 

M omentum : 

(3 H w n 1 _ n 2 U>») (j3 + O + nn M (3H w + n*) C 2 _ ^(3^+**) 3 


= + ^ )2 n(p+n’) 


(2.44) 
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Energy : 

(p+n 1 ) 2 { C (SH^+n*) 2 -^: n' + Sn 2 ♦'} + gnn" (p+ V )(3H w+n *-3) 

= - I^Tt} Fr < 3S w + "♦ ) 2 ( 2 - 45 ) 

The respective equations for no-slip flow situation 
with the induced pressure and its gradient described by a full 
tangent-wedge approximation may also be derived from the 
corresponding equations for slip case. 

Thus, with e^ = eg = 0, equation (2.34) gives the 
momentum equation for no-slip case. 


(3H W n ! - n* O U + YkQT l+b^?+bk]> + n n ' (3H w + mP ) 


Cl- ^ (3H W + i* )C { 2Y M m Cf l+kV+bkj - 


;x7. 


YM 


- — -rr-w 
f 1+b^k 


= Igi M 2 (3H„ + n+ )° i 


2 1 


(2.46) 


The energy equation is obtained from equation (2.43) as 




{ C (3‘H^ + nip ) 2 - 9H W C ri* + 3n ip 1 } { 1 + YkQf 1+b k + bk H * 


t t - t t 

) C 3(H W -1) + ntpC { 2 y M to CT 1+b k + bk 3 


+ nn ( 3H W + n>P 


YM. 


v l+bnc 


> = - 


(T ' 1)M ” * (Si, + T '+) 2 ( 2 - 47 ) 


P.r 


CHAPTER hi 


METHOD QE SOLUTION 

In this chapter we present the solutions to the two 
ordinary second order equations (2.34) and (2.43) with slip and 
temperature -jump boundary conditions. Solutions will also he 
given for equations (2.46) and (2.47) without slip at the 
boundary. For both slip and no-slip cases integration of 
the complete momentum and energy equations is not possible 
due to the lack of initial conditions. Hence an analytical 
solution is provided by solving the strong-interaction equations 
by a series-solution procedure. This provides the initial 
conditions for the numerical integration of the governing 
equations when the full-tangent wedge approximation is employed, 
for obtaining the solution from strong to weak interaction 
regimes. The no slip equations (2.46) and (2.47) are singular 
in nature at the leading edge (S = 0), hence the integration of 
these equations is started from a point slightly away from 
leading edge. However, in slip case the integration may be 
started from the leading edge itself. 

§3.1 analytical sohjtion tor slip case 

The momentum equation (2.26) and energy equation (2.42) 
for the strong— interaction region is now solved by a series- 
solution procedure. 
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We express n and ip as : 

oo 

Ti = 5° 2 a e n<1 

n=0 

OO 

ip = C P 2 b S ns 

n=0 n 


(3.1) 


Therefore, 


i 

ri 


5 


o—l 


2 

n=0 


n 




ri 


€ 


o—2 


2 

n=0 


a.. 


n 


'P 


! 



2 

n=0 


b. 


n 


(a + nq) C nq - 

(a + nq) (a + nq - 1) £ nq - 
(p + ns) ? ns 


(3.2) 


Next, the expressions given by (3.1) and (3.2) are 
substituted in the equations (2.26) and (2.42). Collecting 
the terms of like powers of 5 , we notice that for a compatible 
series and for slip to be a dominant phenomenon at the leading 
edge 


a = 1, q = 1 

and 

p = 0, s = 1 


(3.3) 


Thus , 

2 

n = a Q £ + a-^ 5 + a g £ 

2 

^ = b Q + b 1 5 + bg C 


(3.4) 



32 


Therefore, 


n lg=o " a o 5 n lg=o ~ 2a l 


and 


'5=0 


^ U =0 13 


The coefficients a Q and b Q can he found directly from 
the equations (2.26) and (2.42). we also notice that at % = 0, 
n = 0. 


The momentum equation (2.26) at 5 = 0 yields 


. v2 


. (P+n 1 ) 


22 


[ y(y +1) e u n' - 3(Y-1) H w ] 

(3.S) 


where in expression (3.5) ri' is evaluated at the leading edge. 
The energy equation (2.42) at 5 = 0, in conjunction with (3.5), 
gives 


a. 




2 e ll \ 4 (1 

+ 3 e 2g Pr " 3 1 



2e 


11 


3 e 82 Pr ' 


(3.6) 

Expressions (3.5) and (3.6) jointly give a Q and h o 
as a function of Y , e u , e 22 , Pr and H w ; all of these are 
already specified. 

a x and b 1 are found by equating the coefficients of K 
in the momentum and energy equations, lhus the momentum and 
energy e quat ions (2.26) and (2,42) respectively give 
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a- 


and 


To., = 


a l a 4 “ a 2 a 3 
a l " a 3 


a 4 - a 2 


(3.7) 


where 


a l = r 


1 a i “ a 3 


3Spp f p 

• C 2 g - (Y+l) a e f 2 ] 


a 0 e f " 


a 2 = r 


a 3 = " F" 1 C H 


fc o 11 

- 3I > +1 > a o f2 H 

2 

jj 

Pr 


{ e 2 l(; 3(Y + l)e 22 f 2 + + 4e gg h> 


a. 


a gf 3e i:L g 2 

~ { 4b 0 (*“ pF” + 2h ) +6a o e ll f ^ Y+1 ^ 

Cs (g-f\) - Hi 


(3.8) 


In the above expressions, the various symbols have the following 
meaning : 


r = g {f 2 [6e 11 (f + a 0 )(Y + l) + 24(r+l)e 11 V 24S w^ 


4 r T 


6g[I l+ a 0 e n (7 + 1 )Zl }+ 12a 0 e l a H 1 


s = 18 (y +1) f e\ ± £f 2 H w (3-f-4a o )+g c I] + 96f g h 

I ' 

+ (y+ 1) e 1 < g (1 “ pj 1 ) 


w 


H 


w 


(3.9) 


f = (P + a Q ) 
g = + a o^ + e 22 b o 
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e 2 

11 = - fr> e ll (a o + p) ‘ 


(3.9) 


Thus, the equation (3 .4) supplemented by equations , (3.5) 
through (3.9) gives the leading-edge solution. 

The analytical solutions developed are compared with 
those by Kumar and Jain £24] for the adiabatic slip case. 

For adiabatic surface with ^ = 0 and H w = 1, relation 


(3.5) gives 


M 1 

a o = 11 U=0 = lY+Tfe 


(3.10) 


11 


b o = 0 


and, from (3.7), 


^ = - qfr )2 [4i + 2e u p ] /64 *u- + 811 


P) 


(3.11) 


» 1 = 0 


Kumar and Jain 8 ot tl10 5Efflie TOlueB ° f a ° ^ &1 

in their analysis as those given "by (3.10) and (3.11), which 
serves as a checK on the correctness of the analytical scheme 

deve loped. 

§3.2 MAbYTICAl SO HJ'flOM FOR OASE^ 

a, n • > lfW the basic equations for no slip 
Unlike for slip flow, rne ousio a 

. „ 0 a+ 4-te leading edge of the wedge, 
flow are singular in nature at the xeaan g 

+ c-ome of the initial information 
For the slip case one can extra 
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from the equations itself and hence the series solution 
procedure is comparatively easier. However, for no-slip case 
no such information is available due to the singular nature 
of the equations, which further complicates the already tedious 
job. The equation which we solve here are the momentum 
equation (2.44) and the energy equation (2.45). 


As for slip case, we use series expansion for n and t 
as given by the expression (3.1). 

Substituting for n, n’ , n" , I and 4>' from (3.1) and 
(3.2) in equations (2.44) and (2.45). We find that, for a 
compatible series, 


a 




and 



Therefore, 

n = a 0 s 3/4 + a-L E + a 2 E 5/4 

, —3/4 ■. r— 1/2 v, ^ -1/4 

* = b 0 5 + b l 5 + b 2 -J 


( 3 . 12 ) 


where a Q and b Q in (3.12) are determined from the equations 
obtainined by equating the constant terms in the expanded 
momentum and energy equations. 
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/9 Zt-JLn v,2 , /54 Y-l £ 9 \ „5 , /»x y— t ay £ s - 

— ) a 0 b o + <53 — H W + 53> a o b o + <53 — H w + 34 V a 


81 y-l "2 . 2V - ^4 

o 


- vf?+TJ C a o b o + ^ a o b 0 + = 0 


2 , 2 




(3.13) 


,9 =, 6 k 3 \ ^ /54 w 27v o 5 >,2 
<53 a o V + <S3 H w “ 53 > a o b o 


Itt (9H +11H 2 ) a 4 b + 
64 v w w y o o 


T 54 (E 1)(E 2 ) X (2Y ~ 1)b ° 

L^l (\'“ 1 )( jrI w ) + Y (Y + l)pr 


a 


o 


+ W & i j fe < 6 \ a o b o + ^ a o b o + <> 


0 


(3.14) 


Equations (3.13) and (3.14) are solved for a Q and b Q by the 
Newt on-Rap hs on method. To start the Newt on-Rap hs on iterations 
the range of variation of a Q and b Q are found for the two 
extreme cases of H w = 1 and H w = 0, and is given below : 


H 


w 

0 

1 


o o 

0.87422 0.983014 

1.0813928 0.0 


Once a Q and b Q are found a^ and b^ are determined by 
equating the coefficients gV 4 the expanded momentum and 
energy equations : 

A o B 2 " A 2 B o 
ai = W- A o 


(3.15) 
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where 


k. 


A-t = 


6 aP (|x a b _ H ) - -4 =™y a (6S + 2a b ) 

Y o '-64 oo w' y(y+1) o'- w o o' 


8. 


55 { a o b <f « -S=i a 0 b 0+ 243 ^ ^~± + 57]+ 270 *=i H w+ 225) 


- -2~tp b (6H +2a b ) 
y(y+T) o v w o o' 


Ao — 


w ¥ 0 < ¥ a oV a oV®w ¥ + *w«S ¥ + 12)} 


B o = 55 Ca 0 b o (108H w+ 27a o b o -27) + H w (99H w -63)l 


+ ff¥} 1 4 a o¥ 1 S Vo b o +5 w [ lla o(WoK]> 


a 3 b 
o o 


2, 2 


B 1 = ~T4~ C^V a o b o^ 504H w-- 225 ) + ( 825 V 657 ) H wI] 


+ 

4 


5|§ a 2 H 2 (H -1) + £ 3a 2 b 4 +UH 2 b +12H ab 2 ] 

64 o w v w ' y(y+1) u oo v/ o woo-* 


a 


B 2 = 55 P b o <% b oC^ b ^ 6 ( 3 - 6 V’>oJ + 6V9-UV 

+ ( 108 ^, - 54 )b 0 ) + |Z a 2 p S,, (H w _ 1 ) 

( 3 . 16 ) 

This completes the leading edge solution for no slip 


case 


§3.3. NUMERICAL INTEGRATION SCHEME 


After developing the analytical solutions for leading 
edge we now solve equations (2.34) and (2.43) for slip case and 
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(2,46) and (2.47) for no-slip case, with the full-tangent wedge 
approximation to get solutions which are valid from strong to 
weak interaction regions covering the transition region. 

Integration is started from the leading edge ( 5 = 0 ) 
for the slip case and from a point slightly away from the leading 
edge (C = 0.10) for the no-slip case, A fourth-order Runge- 
Kutta algorithm is used to start the integration. After the 
fourth step, integration is proceeded by employing 'Milne’s 
Predictor' followed by 'Hamming's Corrector' formulae. The 
predicted value is modified by assuming that the local 
truncation error on successive intervals does not change 
appreciably. A provision for the correction of truncation 
error is also given. These modifications and corrections enhance 
the convergence. Step-size at the start is taken as 0.025 and is 
increased by tenfold after every eighty steps. Also after 
every eighty steps, the kunge-Kutta method is used again to 
proceed through the next four steps and thereafter is followed 
by the predictor-corrector method. This procedure avoids the 
chances of the errors getting magnified and also enables a 
faster convergence. 

The various formulae used in the numerical integration 

are : 

(a) Pourth order Melne' s predictor formula : 

y i+l,0 = y i-3 + I h ( 2f i " f i-l + 2f i-2) 


(3.17) 
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(Id) Modification for the predicted value is obtained from j 

1 12 

n + i )0 = n + i,o + m (yi )k - n, 0 ) ( s - 10 ) 

(c) Hamming's Corrector equation 

y&i = 5 { 9yi-y^ sll Cf(x 1+1 ,ii + i +2f i- f i-i3> < 3 - 19 ) 

^ — Oy-^s' :9 5***9( k— 1 ) • 

(d) Correction for truncation error is given by s 

^i+1 ~ ^i+l,k “ 12 1 ^i+l,k “ ^i+1,0^’ (3.-^0) 

where ! f' is the gradient of y in the direction in which the 
step 'h* is taken. Here * i* is the step count and ' j' is the 
iteration count. 

The solutions were obtained for various values of the 
free stream Mach number, semi-wedge angle, wall enthalpy ratios 
and Prandtl number. 

Computations have been carried out on the IBM. 7044/1401 
computer at the I.I.T. Computer Centre. It is found that one 
set of computations requires a computer time in the range 1,5 
to 2 minutes. On the other hand, a finite difference technique 
j^26^ requires a computer time not less than 30 minutes for 
one set of computations. Thus, the present method saves the 
computational time considerably. 



CHIP TER IV 


RESULTS AWL DISCUSSIONS 

In this chapter, the various results obtained in the 
present investigations and their comparison with the previous 
theoretical and experimental findings are discussed in detail. 
Briefly, the conclusions drawn therefrom are presented at the 
end of this chapter. 

Pig . 2 shows the variation of the boundary layer 
thickness along a cold surface for semi-wedge angles of 0°, 

2° and 5°. It is seen that the boundary layer thickness holds 
an inverse relation with the wedge angle, which is true because 
a higher wedge angle causes an increased inviscid compression 
which, in turn, will have a tendency to keep the boundary 
layer thin. It is also seen that the boundary layer thickness 
is larger for Ivi^ = io than for If M = 2o for the same values of 
U and the semi-wedge angle. This behaviour is also attributed 
to the higher inviscid pressures at high Mach numbers which 
overcome the opposing effect of rarefaction at high speeds. 

Our results are in quite good agreement with those of 
Ghatt op adhyay 27^j for semi-wedge angles of 0° and 2° and 

for Mach numbers 10 and 2o for no-slip case. 

In fig. 3 a plot of 6 for various thermal conditions 
at the wall for both slip and no-slip cases is given. It can 
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"be seen that 5 continues to decrease as the wall becomes cooler. 
This is expected because the density of the fluid in the boundary 
layer is maximum when the wall is fully cold. Again, for the 
same wall-enthalpy ratios, slightly lower results are obtained 
when the effect of rarefaction is incorporated. It is worth 
noticing that the difference between the slip and no-slip 
solutions becomes negligible at highly cooled surface conditions. 

The surface pressure distribution as a function of x 
is shown in fig. 4 for different wall temperatures. We observe 
a reduced pressure level for a lower wall temperature. Similarly 
slip reduces the pressure level markedly towards the leading 
edge and much less at the downstream stations. This type of 
behaviour is due to the fact that, in a hypersonic stream, the 
effective thickness of the body decreases with decrease in wall 
temperature or slip. In fig. 4, present results are given for 
three wall conditions H = 0.5, 0.15 and 0.065 with M = 20, 

Pr = 0.V5 and (3 = 0° for slip and no-slip cases. For H = 0.065, 
present results are compared with the experimental data of 
Vidal and Wittliff £28^ and for = 0*15 a comparison is 
made with the data of Vas and McDougall on. A reasonable 
agreement is observed in these comparisons. The theoretical 
predictions of Gheng £ 13 ^ for strong and weak interaction 
zones agree quite well with the present curves in the relevant 
regions. However, the present results are slightly lower than 
Menon’ s OH finite difference solutions.’ Pig. 5 shows the 
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variation of wall pressure along the surface of a semi-cooled, 
body for various wedge angles. It is seen that the pressure 
level increases with increase in wedge angle and Mach number. 

The results compare favourably with those of Ohattopadhyay 

- 3 - 5/2 

The variation of 0^/2 x with rarefaction 

parameter is shown in fig. 6 for four different cases of wall 
cooling, namely, = 1.0, 0.5, 0.20 and 0.065. The skin 
friction is seen to decrease as the wall is progressively cooled, 


due to the reduction in the coefficient of viscosity. For the 


same wall conditions a lesser skin friction is noticed for the 


slip case. This is obviously due to the lower velocity gradient 
at the surface. Results are shown upto ¥ 2 where the flow is 

molecular in nature and the present theory does not hold good. 

At ¥ 2 l 0.4 the present theory deviates considerably from the 
experimental data of ¥idal and Barts £ 20 "2 • Theoretical results 
of Shorenstein and Probstein [^22^, Rudman and Rubin [19] and 
Suresh Menon £ 26 ] are slightly higher than those obtained in 
the present work even for ¥ < 0.4 as shown in fig.- 6. Pig. 7 
shows the effect of semi-wedge angle on skin-friction coefficient. 


The skin-friction coefficient increases with increase in wedge 


angle for the obvious reason of larger velocity gradient at the 

wall. The no-slip solutions give constant values of M^0^/2x / 

for ¥ around 0.10. Por the slip case on the other hand, the 
3 _3/2 

quantity M m G^/2x decreases markedly towards the leading 
edge . 
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3 ■ ,- 3/2 

Variation of the quantity M m St/x as a function of V 

for various values of H^. is shown in fig, 8. Stanton number is 

observed to be increasing with increase in wall enthalpy ratios. 

The present results compare quite favourably with the experimental 

data of Vidal and Bartz £ 20 3 . Bor V above 0.4, a comparison 

with the results of Rudman and Rubin £ 19 is also given and 

there is good agreement. Although the theoretical predictions 

of Shorenstein and Probstein £22^j do not agree with the present 

work above V = 0.4, they show a reasonable agreement with present 

theory in the downstream limit. Me non £ 26 also gives similar 

results for the merged layer and downstream regimes. In fig. 9, 

curves (1), (2) and (3) show the effect of wedge angle on the 

3 ,- 3/2 

variation of the quantity St/x as a function of the rare- 
faction parameter. The Stanton number increases with increase 
in the wedge-angle, which is also indicative of the increased 
heat transfer. The curve of St/x 3 ^ for = 20 falls below 
that for IVI = lo for the same values of H , Pr and 8. A 
comparison between curves (1) and (7) shows that for constant 
values of E^, and {3, the Stanton number decreases with 
increase in Prandtl number. 

The variation of slip velocity with rarefaction parameter 
on a flat plate for different wall enthalpy ratios is shown in 
fig. 10. Results are shown up to V 3 although these extreme 
regions fall beyond the scope of the continuum analysis. Slip 
is shown to be clearly dependent on wall temperature, tending 
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to zero at the extreme cold wall case. Towards the leading 
edge, the slip velocity becomes a considerable fraction of the 
free stream velocity. The present results compare reasonably 
well with the experimental findings of Vidal and Barts £2o[] 
and the theoretical predictions of Shorenstein £22^ and 
Me non. £ 26 ^ upto V 1.0.3 and thereafter shows a comparatively 
slow rate of increase. Pig. 11 shows the effect of semi-wedge 
angle and Mach number on slip velocity. Slip is seen to decrease 
with increase in wedge angle at constant wall cooling. A 
comparison between curves (1) and (4) clearly shows the 
dependence of slip velocity on Mach number. The slip velocity 
over an insulated flat plate is also shown in the same figure 
and an excellent comparison is noted with the integral analysis 
of Kumar and Jain |_24[] . It is interesting to note that the 
analysis shows the effect of slip even for V as low as 0.01, 
v/hich is a region well within the scope of the strong interaction 
theory. This figure also gives a first order estimation of 
(c. which is approximately equal to ~ . Por the strong 

Q U oo 

interaction regime, this quantity is much less than unity. 
Therefore the approximation (2.17a) is justifiable. 

Pig. 12 shows the variation of enthalpy jump as a 
function of rarefaction parameter for different thermal 
conditions on a flat plate. It can be seen that the enthalpy 
jump increases as the gas becomes more rarefied. Also a heated 
surface indicates an increased enthalpy jump in comparison to a 
cooled surface due to the greater thickness of the Knudsen 
layer at the plate surface in the former case. In the range 
0.055 < V 0.4034, the results obtained in this work show 
reasonable agreement with the experimental data of Vidal and 
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Bartz £ 20 ^ . Theoretical investigations of Shorenstein 
and Probstein £22]] and Menon £26£j, however, show a 
slightly lower enthalpy jump . 

Pig. 13 shows the temperature jump as a function of 

X for different levels of wall cooling. Towards the leading 

edge, cooling the wall up to = 0.6 gives negative temperature 

jump. Purther downstream, the gas rapidly adjusts to the wall 

temperature. Por wall conditions like H = 0.2 and H = 0.4 

w w 

the temperature jump is positive even close to the leading 
edge. Thus the temperature jump can he positive or negative, 
depending upon the wall enthalpy variation. This phenomenon 
was studied in detail by Metcalf £ 30 J . 

The enthalpy profiles at various plate locations for 
a semi-cooled wedge are shown in fig. 14. It is seen that 
the profile becomes fuller as the leading edge is progressively 
approached. Pig. 15 shows the effect of wall cooling and 
semi-wedge angle on the distribution of enthalpy across the 
boundary layer at x = 0.5. Curves (1), (2) and (3) indicate 
that the enthalpy profiles become fuller as the wall is heated. 
Pig. 16 gives the velocity variation across the boundary layer. 
It is seen that the velocity at the wall decreases with the 
decrease in wall temperature and strong interaction parameter, 
and is virtually zero for the fully cooled surface. 



CONCLUSION 


In light of the detailed comparativ e studies carried 
out with the experimental as well as theoretical investiga- 
tions 3 it is seen that the approximate integral method is 
quite powerful to handle the hyper sonic viscous-inviscid 
interaction problem with a reasonably good accuracy . The 
method is relatively simple to use and requires a lesser 
computational time compared with other techniques . The 

formulation is so generalised that it includes all the possible 
flow situations and the solutions appear as a function of 
Mach number 3 Reynolds number 3 Prandtl number 3 Knudsen number 3 
wall temperature and semiwedge angle. 

Analysis shows that slip exists even in the strong inte- 
raction region and tends asymptotically to about 90% of the 
free stream velocity towards the leading edge region. The 
rarefaction effects are reduced with increase in wedge angle 
and wall-cooling . For the fully cooled surface 3 slip is 
apparently absent. The slip solutions give lower values 
of surface pressure 3 skin friction , boundary layer thickness 3 
heat transfer etc . 3 compared to no slip solutions. Another 
noticeable observation is that the temperature jump at the 
surface can be positive or negative depending upon the value 
of wall enthalpy ratio. 
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Thus we can say that the slip and temperature jump 
boundary conditions at the surface are important to explain 
the experimentally observed plateaus of constant pressure, 
skin friction , heat transfer etc. and in obtaining results 
closer to the experimental observations . 
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APPENDIX 


§A-1 THE SLIP PARAMETERS 


Prom the kinetic theory (31) we may write : 


X m = 1.256 Y 


1/2 


or 


1 w 


A w 

IT 


= 1.256 Y 


1/2 


u /a p 
H w' w w 


P„ u 


% p 


w 


= 1.256 y 


1/2 


U oo 

a oo a -^ P W 


= 1.256 Y ^ "J-oo V 7JT 


W Pw 


x.e 


A w 

r ~ 


i.256 yy 

^7p7 


Further, 


T H 

w w 

I “ " 1 2 

00 H oo - Z ^ 


= \ (! + ^ M «)* 

2 

For hypersonic flows with M >> 1, one can 


(^i) M* 


Substituting (A-2) in (A-l) gives 


(A-l) 


write 

(A-2 ) 



A2 


X w 1.256 M 


2 J Y(Y-i; 


(P/Poo) 


(A-3) 


When pressure is described by a tangent-wedge approximation 
valid for strong-interaction region only as given by(2.22), one 
obtains from (A-3) : 


Xw 1.256 

Tj ~ Y+1 




r~o — asw ' 
( P + cE> 


(A-4) 


Equations (A— 4) enables us to define the two slip 


parameters as 

Cl X, 


= (P W) 2 = ^ Cl [2^i w 


c 9 x 
2 w 


. s2 1,256 


CP + n ) = “5 0 2 8^H w 


(A-5) 


We also define two other slip parameters from (A-3), 
such that it doesn’t impose any restriction on how the pressure 
is being described : 


e l = 


0-1 X 1 
1 w 1 


- = 1.256 O-i £ y Hi 

CO X •— C W •—* 


Lull H 1 V2 


0 o X -i 
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wherein (A-6 ) 


» = «C/(p/p j . 
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§ A.2 . "VELOCITY AND ENTHALPY PROFILES 

Eor computational purposes the velocity and enthalpy 
profiles given by (2.12) and (2.16) respectively are non- 
dimens ionalised using the characteristic length ’ L* as given 
"below: 


Equation (2.12) gives 


u 

C- X + y 

1 w J 

U oo ~ 

c i v + ^ 


^1 *W / T 

— E — + y/L 


0-1 X C 

1 w . 6 

I + x 

U 

6 1 (U + 4 

*oo " 

o) + n 


(A-7) 


wherein- (A-7),$ = y/1. 

Similarly the expression (2.16) gives the n on- dimens ionalise-d 
enthalpy profile as 


2 , 2 


H w = (H w + e g oji(i) + *<► + £ U-H w ) - i|»(e g u + n)I] /n 


wherein (A-7) and (A-8), $ varies from zero to n 


(A-8) 


The slip velocity, u^ is computed form (A-7) with <)> = 0 : 


u. 


e 2 o) 


u co e i “ + n 


(A-9 ) 
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§ A.3. SKIN FRICTION GQBHi'IGIBNl 

By definition the shear stress at the wall is 


T W ' ?,M'\ 


l w ^ 6y ; w * 


Substitution for (~) w from (2.12) leads to 


u u. 

H \Y oo 


w °1 »w + 5 ■ 


Now, the skin friction coefficient may be defined as 
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0 f " 2 p~ FBI *„ + 6) - 0, x_ + 6 


)s\ /V 

1 w 


or 


G* = 


f to + ri 
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5A.4. STRONG DI1BRAC1IQN PARAMATER 


The strong-interaction parameter is defined as 


X = M oo v — • 
X 


further, using the linear temperature-viscosity law 


t 2 

— = c in the following definition of 5 
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15 


We obtain, 
Re 


x 

T 


Re 


w 


s w d + ¥ m ~> c 


or 


M 


5 = 


(“® 1 o/Re x ) 8 S^l + Ifi «£) 


MJ 


XT /I Y— l-nirSv —2 
H W ^ + “2" M J X 


This may be rearranged to give : 
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§1.5. HEAT TRANSFER AND STANTON NUMBER 

By definition, the heat transfer at the body surface 
is given by 

% = If). • 

Also, from the definition of H : 

/ & T \ __ 1 / dH \ 1 / du. \ 

v '9y'w ~ TT~ '•dy'w “ S’ ^3y 'w .* 

P 

2 

Evaluating (§^“) w an(i (|~) w from the expressions (2,12) 
and (2.16) for the velocity and total enthalpy profiles, 
respectively, we get 



A6 


% 




k H„ 


: ) r b - 

'w [_ 


20 1 




or q, 
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(A- 12) 


Stanton number, by definition, is 


St = - 


or 


S.'s 


P oo U oo ( B o “ 


Poo U oo H m (1 - H w ) 


Substituting from (A- 12) gives t 


St = 


s r ( 1 - V 




(A- 13) 


All these results for no-slip cases can be computed 
by equating the slip parameters to zero. 



